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rS^ ' Abstract. We propose an implementation of a valley selective electronic Veselago 

j/3 ■ lens in bilayer graphene. We demonstrate that in the presence of an appropriately- 

oriented potential step, low-energy electrons radiating from a point source can be 
re-focused coherently within the same band. The phenomena is due to the trigonal 
warping of the band structure that leads to a negative refraction index. We show that 
the interference pattern can be controlled by an external mechanical strain. 
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1. Introduction 

Since its first experimental isolation [HE] graphene has had a profound impact on 
the solid state community. The past year has seen reinvigorated experimental and 
theoretical investigations of bilayer graphene, as sufficiently-clean samples have been 
fabricated to probe the low-energy spectrum [3]. A key conclusion of these investigations 
is that mechanical distortion as well as electron-electron interaction can lead to profound 
changes in the topology of the band structure and to symmetry breaking in the 
electronic system [3H6]. In recent years, new ideas in electron optics have utilized the 
exotic electronic structure of graphene and other novel materials such as topological 
insulators. Electron lenses utilizing negative refraction [7HTU] and valley-polarized 
electron beam splitters [TTHTi] among others have been proposed to manipulate electron 
beams in mesoscopic devices. Apart from the manipulation of the valley degree of 
freedom, graphene devices have already been proposed for spin polarized electron beam 
engineering using magnets [15] , and for spatial separation of electrons and holes using 
normal-superconducting-normal junctions [16J. 

Unlike 2DEGS (two-dimensional electron gas) in more conventional semiconduc- 
tors [I71IT8], which require a magnetic field to provide the focussing of electrons and 
imaging of the interference pattern of injected electrons is hindered by the fact that 
the 2DEG is often buried inside a heterostructure [15] . graphene-based devices can be 
self-focussing and their interference patterns are exposed on the surface. In the present 
manuscript we propose an electron-optical device, whose valley- dependent interference 
pattern can be controlled by strain. By analogy with conventional 2DEGS, where the 
control of focussing by an external magnetic field leads to sensitive magnetic field sensors 
even at room temperature |20j, our graphene-based device forms a basis for a sensitive 
detector of strain. 

2. Electron Focusing in Novel Materials 

Focusing of electrons in graphene by a planar pn junction has first been proposed by 
Cheianov and coworkers [7J. In their proposal, particles are focused due to the fact that 
in the valence band, the group velocity points in the opposite direction to the wavevector. 
Hence electrons emitted from a point source on one side of a pn junction converge on the 
other side at a focal point, thereby realizing a Veselago lens [2T]|22] . Circular pn junctions 
have also been investigated recently in single layer [8] and bilayer [23] graphene, where 
apart from the focal point, a hierarchy of caustics was discovered and described using 
the semiclassical approach of catastrophe optics [21]. Later it was pointed out that 
warping of the band structure [25J can also induce a focusing effect even when particles 
on the two sides of the junction are in the same band [9j. This was demonstrated with 
the hexagonally warped surface states of the three dimensional topological insulator 
Bi 2 Te 3 [26J. In this work we propose a flat interface setup (shown in Fig. [I]) that 
can focus electrons in bilayer graphene due to the trigonal warping [25] of the band 
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Figure 1. Single band focusing in bilayer graphene with a potential step. A particle 
emitted from a source disperses or focuses depending on which valley it came from. 
Measuring the maxima of the transmitted electron density with an STM tip can be used 
to indirectly gain information about the low energy topology of the band structure. 



structure. In this section we review some aspects of electron optics with regards to 
focusing particles by a flat interface. 

For a particle incident on the barrier with wavevector k« and group velocity 
v i = (•9k-^(k)//i)| k=k ., the refractive index n at the barrier is defined by the states that 
conserve the momentum parallel to the interface and propagate away from the interface. 
Let k t/ y and v t / r denote the wavevector and group velocity of the transmitted/reflected 
particles respectively. If the velocities Vj/ 4 make an angle <pi/ t with the normal of the 
barrier then the refractive index is 

«) 



sin 



n 



(1) 



smWt) 
In general, n depends on the direction of k;. 

The intraband focusing effect of a planar interface is closely related to the concave 
geometry of the £'(k)-energy contours of the dispersion relation. It can be shown that 
electrons radiating from a point source into a small region in k-space will be focused 
on the other side of the junction if and only if the curvatures of the energy contours on 
the two sides at the corresponding k-wavevectors differ in sign. If this holds around the 
k-vectors aligned with the optical axis, that is for rays that are nearly perpendicular 
to the interface, then the different signs of the curvatures result in a negative refractive 
index, which is however not a necessary condition for focusing if we move away from the 
optical axis. In particular, if the _E(k)-contours are locally symmetric to the direction of 
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Figure 2. Real space trajectory and fc-space directions of a single particle incident on 
a potential step in the valley K' in a bilayer graphene junction at u> = and a = 0, 
i.e. where the (x, y) and (x',y') systems coincide. Angles op denote the direction of 
the wavevectors, while angles <f> stand for the direction of the group velocities. Note 
that in this anisotropic system the group velocity and the wavevector of a given state 
are in general not parallel. A particle emitted from the point source S at — 2.07/im 
on the x-axis, hits the boundary at the y-axis with an angle of incidence <f>i and along 
with electrons on neighbouring trajectories, gets diverted towards a cusp type caustic, 
the maximum of which we call the focus F, for the sake of simplicity. The electrons' 
energy is 80 meV on the left and 10 meV on the right side of the junction. Blue 
arrows represent the incident particle, green arrows the transmitted, and red arrows 
the reflected ones. The wavevectors are denoted by dashed arrows and the classical 
trajectory of the particle by solid lines. Solid lines depict the electrons' dispersion 
curves on both sides of the junction. Note that the curvature has different sign on the 
two sides for particles close to the optical axis. 



the optical axis, then a simple, direct relation can be formulated between the refractive 
index for these rays and the curvatures on the two sides [9]: 

_ c(kj, Ej) 

c(k t , Et y {2) 

where E{ is the chemical potential on the incoming side of the junction and E t is the 
chemical potential on the other side of the junction. It also follows that around the 
optical axis, the negative sign of n is needed for the focusing phenomenon. In Fig. [2] 
we show the trajectory of a single particle hitting a planar potential step in a bilayer 
graphene sample. (The barrier coincides with the y'-axis.) A classical trajectory can be 
parametrized by ipi which determines all the remaining angles. 

If a point source is situated at (— x s ,0), then on the other side of the interface, a 
generic trajectory has the form 

y = x s tan (^ (</?*)) + x tan (0 t (<#)) . (3) 

Some of the refracted rays touch each other at certain points. The sets of such points 
are called caustics, which are envelopes of a family of rays at which the density of rays 
is singular. In the classical theory of geometrical optics, this singularity is unphysical, 
however, the results of the semiclassical catastrophe optics [27] agrees well with the 
quantum mechanical treatment, which predicts a finite particle density with local 
maxima [21]. The focal point of particles corresponds to the global maximum of a 
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cusp type caustic, which evolves into two fold type caustics. The coordinates of these 
curves can be obtained by differential geometry: 

dnV = 0- (4) 

In order to describe the system quantum mechanically, one has to solve the 
Schrodinger equation on both sides of the junction and match the solutions at the 
boundary. The boundary conditions can be derived from the Schrodinger equation by 
applying it for an infmitesimally narrow stripe containing the junction. In the case of the 
model Hamiltonian used in this work, these conditions require the wavefunction to be 
smooth and the same for its first derivative with respect to the direction perpendicular to 
the junction. This can be satisfied with 2 planewave-like solutions of the homogeneous 
Schrodinger equation on either side of the junction. For a given k y /, which must be 
conserved during the scattering process, there always exist 1 or 2 propagating waves and 
1 or decaying wave, giving always exactly 2 solutions on either sides needed to satisfy 
the boundary conditions. By solving this system of equations, we obtain the coefficients 
for the reflected and refracted rays matched to the incident ray. Since all the calculations 
can be extended in a straightforward manner for several rays, for the sake of simplicity, 
wherever applicable, we will assume that we only have one refracted wave propagating 
to the right direction. To calculate the density of particles, the contributions of all the 
refracted rays with the corresponding transmission amplitudes were summed up in the 
far field approximation. By the summation of the rays, a homogeneous distribution of 
the source electrons in ipi is assumed. 

3. Focusing in Bilayer Graphene 

The central idea of our work is to use bilayer graphene with a potential step for 
focusing electrons, where particles remain in the same band. As we will demonstrate, 
this idea retains all the neat features of earlier setups [3 [9], including focusing and 
high transmission, and it even incorporates the idea of valley- dependent electron beam 
manipulation [TTUT21I26]. 

We start with the Hamiltonian of the system under consideration 

H _( V{r) -±-p 2 _ + & 3 p + + w\ 

\-£pI + &3P- + w* V(r) J' {) 

where p± = p x ±p y are the momentum operators and £ = 1 for valley K and £ = — 1 
for K'. The material parameters of bilayer graphene m ~ 0.035m e and t>3 ~ 10 5 m/s 
are defined among others in Ref. [28 J. For very low energies (< 1 meV), the complex 
parameter w describes symmetry breaking due to the electron-electron interaction [I], 
but it also describes the external mechanical strain [5] for higher energies as well. The 
relation between w and the mechanical distortion of the lattice can be formulated as 

w = ri(8-6')exp(-2i9), (6) 
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where 77 is a real constant describing mechanical properties of the lattice, 5 > 5' are 
the principal values of the strain tensor \{didj + djdi), and 9 is the angle between the 
principal axes and coordinate axes (x,y). In the strain tensor, % and j are the (x,y) 
coordinate- indexes and (d x , d y ) denote displacements from the equilibrium positions. In 
the expression (jHJ) we neglect the interlayer shear shift. The potential V(r) describes a 
potential step that is sufficiently smooth on the atomic scale so that it does not introduce 
inter- valley scattering, and we can also neglect a term proportional to the gradient of the 
potential that arises in the Schrieffer- Wolff transformation in the derivation of (jSJ) [25JT29] . 
We assume V(r) to have the form 

V{r)=V 1 + V 2 e(x-ta,n(a)y), (7) 

i.e. the potential on the two sides of the step is constant V\ and V2 and the barrier 
makes an angle a with the crystallographic TK direction (A^-axis). It is convenient to 
introduce a coordinate system (x',y') where y' is defined by the planar potential step, 
and the source is sitting on the axis x', which means that x and x' and also y and y' make 
an angle denoted by a. We do not take into account a gap term in the Hamiltonian 
that accounts for the electrostatic asymmetry of the two graphene layers, since this 
would not alter the main message of the present work. We also neglect second nearest 
neighbour hopping in the individual graphene planes that has negligible effect on the 
results presented here. In the absence of the potential V(r) the spectrum of (jSJ) has been 
studied in [51125] . For low energies, depending on the value of w, the spectrum consists of 
four or two massless Dirac cones, whereas at higher energies, the band structure merges 
to form only a single pocket per valley. A typical energy scale of this model is the 
Lifshitz energy El = mv%/2 ~ 1 meV, which defines the Lifshitz transition of the band 
structure's topology from multiple pockets to a single pocket in the absence of strain 
(w — 0). Since the Lifshitz energy is quite small, it was difficult to gain experimental 
data from energies below it until recently [3] . 

As shown in the previous section, focusing is controlled by the curvature of the 
band structure. If on the two sides of the potential step there exist regions in the 
A;-space on the given energies with overlapping projections to the dimension parallel 
with the junction, and for these corresponding fc-regions the curvature of the dispersion 
relation has different sign, then the system can focus electrons. In Fig. [3] we present 
the curvature c of the E(k) dispersion curve of the considered model as the function of 
energy along the TK direction for different values of w. It is clear that the curvature 
in the absence of w changes sign roughly at E « 42.1 meV. This energy scale is still 
sufficiently small compared to the direct inter layer hopping [251I2E] 71 ~ 0.4 eV, that is 
it falls in the energy range where the two band model (jSJ is still valid. Since a changing 
w can produce a large change in the energy of zero curvature, one can study the focused 
electrons in bilayer graphene samples with a planar potential step, and yield information 
about the value of w, even in samples that are only moderately gated (~ 50 meV). 

Our setup also has the capability of separating electrons in different valleys, without 
the need for an unrealistic potential barrier height of around 1 eV, as in previous 
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Figure 3. Curvature c as the function of energy along the direction indicated by 
the dashed line in the inset. Inset shows the energy contours around the K' point 
for w — meV. In the main figure, solid line denotes curvature for w = —6 meV, 
dashed line is for w = meV and dotted line is for w = 6 meV. Black dots denote the 
particles' energy on the left (80 meV) and right (10 meV) side of the junction as in 
Fig. [3 The electrons' energy is chosen to find optimal conditions for focusing, that is 
a negative and slowly varying refractive index n in a possibly wide range of the angle 
of incidence 6. 



graphene based electron optical setups (THE]. This is demonstrated in Fig. HJ where 
we follow the classical trajectories of several particles emitted from a point source, in 
real and in fc-space for the case of undistorted samples, with the interface aligned on 
the ?/-axis i.e. w = meV and a = 0. It is clear from the figure that electrons from the 
two valleys behave differently. Particles in valley K disperse upon passing the barrier, 
while those in valley K' gather in a focus point on the other side. If the barrier is not 
aligned with the crystallographic orientation of the lattice, then in general there will be 
two or three foci each gathering particles from different valleys as shown in Fig. [5] 

One key ingredient that is not captured by the ray tracing analysis is the 
transmission probability of particles reaching the junction. This can be calculated by 
evaluating first the expectation values of the current operator along the x'-axis. Then 
the transmission probability T will be the ratio of the current carried by the refracted 
wave in the x' direction to that of the incoming wave also projected to x' . The reflection 
probability R can be calculated on a similar basis from the reflected wave. Note that 
the evanescent waves do not carry current. 
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Figure 4. Real space trajectories and fc-space directions for both valleys with barrier 
parameters same as in Fig. [5J In figures b) and c), blue arrows show the direction of 
the group velocity of incident particles with wavevectors corresponding to the arrows' 
position on the -E(k) energy contour (grey solid curves). Green arrows are for the 
transmitted, and red ones for the reflected electrons in both valleys. The arrows in 
a) and d) show the real space trajectories of particles emitted from the source (blue) 
and partially transmitted trough (green) or reflected by the junction (red). Note that 
particles in valley K' focus while particles in the valley K diverge. We only show 
trajectories that have a finite probability of transmission at 80 and 10 meV. 
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Figure 5. Total particle densities and theoretical curves of the caustics (blue lines) for 
different values of a. Note that a given focus point only gathers particles from a single 
valley. The position of the point-source is at —2.07 /im on the x'-axis. The electrons' 
energy is 80 meV on the left and 10 meV on the right side of the junction. The equation 
of the sheared cusp caustics can be derived from equations ([3]) and fU) 1HJI30] ■ 
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Figure 6. Quantum mechanical transmission probability as the function of angle of 
incidence in valley K' 1 for different values of w as in Fig. [3] The barrier parameters are 
the same as in Fig. [2] 



In Fig. [6] we show the transmission probability as the function the angle of incidence 
for different values of w in the case of a = 0. Remarkably, the transmission probability 
is high and roughly independent of the angle of incidence fa for particles reaching the 
junction at a wide range of angles of incidence around the x'-axis. This is a beneficial 
result of scattering within the same band. For particles that change band during 
refraction at the junction, the transmission probability is reduced, is very sensitive 
to fa, and vanishes exactly for perpendicular incidence [23], unless gaps are induced by 
external gating [10]. 

For finite a, when the junction is not aligned with the y-axis, the focus point will 
not remain on the x-axis but will continuously move with a. This is depicted in Fig. [5] 
where the total particle density is plotted with contributions from both valleys and for 
different values for a. We also calculated the classical curves of the caustics using (jJJ). 
In the figure, these curves are represented by blue lines, the crossing points of which 
approach neatly the maxima of the particle density as expected. There are two or three 
focus points each gathering particles either from the K, or the K' valley. 
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Figure 7. Calculated density of transmitted particles for different values of w. The 
particles' energy is 80 meV on the left and 10 meV on the right side of the junction. 
The source is at — 2.07/im on the rr'-axis. The sub-figures correspond to a) w = meV; 
b) w = 6i meV; c) w = — 6i meV; d) w = —4 meV; e) w = 2 — 2i meV. 



4. Effects of broken symmetry 

As shown in Ref. [I], the electron-electron interaction in bilayer graphene leads to a 
spontaneous breaking of the threefold rotational symmetry giving rise to a finite value 
of w in the effective low-energy Hamiltonian (J5J). However in this description of the 
electron-electron interaction, the Fermi energy as well as the energy of the excitations 
around it need to be of the same order of magnitude as the Lifshitz energy (~ 1 meV), 
which is not the case discussed in this paper. Also, w originating from the interaction 
can have very different values on the two sides of the junction, which are unknown. 
On the other hand, later, it was argued in Ref. [5] that mechanical distortion leads to 
the same low-energy effective theory, with a constant w on a wide energy range. A 
key consequence of a finite w is that depending on its value, the low energy spectrum 
of the system can undergo a topological phase transition. For small magnitudes of w, 
the spectrum consists of four Dirac cones, while at larger w, there are only two Dirac 
cones in the spectrum. The small value of the Lifshitz energy makes it difficult to infer 
directly the number of Dirac points. Therefore it is essential to seek indirect methods 
to determine the value of w and hence the number of Dirac points. 

As presented in Fig. [3j the curvature of the energy contours depends considerably 
on the value of w. This suggests that the position of the focus point will be similarly 
sensitive. Indeed, this is clear from Fig. [7] where we show particle densities from valley 
K' for different values of w. The most important information one can extract from this 
data is the position of the maxima that can be used to directly infer the value of the 
complex parameter w. In Fig. [BJ we show the semi classical estimation for the position 
of the tip of the cusp caustic as the function of the real and imaginary part of w for 
a = 0. From the expression ([6]), it is apparent that the imaginary phase of w will be 
—26 reflecting a direct relation to the geometrical deformation of the lattice. It is also 
worth mentioning that if 9 is a multiple of tt/2, then w is real, but if mod(#, vr/2) = 7r/4, 
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Figure 8. A semiclassical estimation for the position of the tip of the cusp caustic 
formed by electrons from the valley K' as the function of the parameter w for a = 0. 
The electrons' energy is 80 meV on the left and 10 meV on the right side of the junction. 
The position of the source is at —2.07 /im on the x-axis. Inside the green curve the 
low energy spectrum has 4 Dirac cones, between the red and green curves there are 
two Dirac cones and a local minimum, while outside the red curve the spectrum has 
two Dirac points but no local minimum. The nonzero value of w is attributed to a 
distortion of the lattice that is around or less than 1% on this map [5j. Here we note 
that the e — e interaction renormalizes w, thus rendering it slightly energy dependent 
and changing its relation to the applied strain. This however does not affect our results 
qualitatively. 



then w is purely imaginary, and otherwise it is a general complex number. Interestingly, 
if 5 = 5', then w vanishes. This corresponds to a 'hydrostatic' rescaling of the lattice 
and only affects the values of m and t^. Based on the results of Ref. [5j and [51] . it can 
be show that a distortion of 1% changes m and t>3 just by a few percent and thus the 
effect on the diffraction pattern is negligible. On the other hand, we have shown that 
the position of the foci sensitively depends on the value of w. Thus by experimentally 
investigating the focusing effect, one can measure the external strain present in the 
sample. 

In clean systems, the most important factor that determines the lifetime of 
excitations is electron-phonon scattering. According to recent calculations j32j|33], the 
mean free path for excitations in the considered energy range around 50 meV is of the 
order of microns, which is accessible by current experimental techniques |34j . 
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5. Conclusion 

In this paper, we investigated the anisotropic electron optics of bilayer graphene. We 
demonstrated that a moderate potential step can be used to focus electrons within 
the same band in a valley-selective manner with high transmission probability. We also 
investigated the effects of broken symmetry in the electronic structure due to mechanical 
distortion. The presented results clearly show that the proposed device can be used to 
determine symmetry breaking and extract indirect information about the low-energy 
topology of the band structure of bilayer graphene samples. These results are expected 
to form key ingredients in future device design, where one can envisage using bilayer 
graphene as building blocks for more complex electron-optical devices. 
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